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ABSTRACT 
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I. INTRODUCTION 


Boc OBJECTIVE 


Advanced technology propulsion systems require a control 
system design that is able to achieve the desired performance 
at the greatest possible efficiency. One particular scheme 
that shows great promise in this regard is the optimum integral 
control design program, CONSYN, developed in Ref. 1. 

This study was undertaken to examine the validity and 
effectiveness of applying CONSYN to the FFG-7 ship class gas 


turbine propulsion system model obtained from Ref. 2. 


B. COMPUTER FACILITIES AND SOFTWARE UTILIZED 


The above objective was implemented by simulation on the 
IBM 360/67 computer of the W. R. Church Computer Center at 
the Naval Postgraduate School. The following programming 


software was required: 


e 


Various routines in the International Mathematical and 
Statistical Library (IMSL) 

" The IBM developed Continuous System Modeling Program 
(CSMP) 

The Control Program for Engineering Synthesis with 
Constrained Function Minimization (COPES/CONMIN) de- 
veloped by Garret N. Vanderplaats of the NASA Ames 


Research Center 








The routines developed for CONSYN in Ref. l as modified 


by Ref. 3 and this study, are listed in Appendix A. 





II. PROPULSION PLANT MODEL 


The propulsion system for the OLIVER HAZARD PERRY (FEFG-7) 
class patrol frigates consists of two General Electric 
Corporation LM2500 Marine Gas Turbines, coupled through a 
reduction gear to a controllable reversible pitch (CRP) 
propeller manufactured by the Bird-Johnson Company. Since 
the gas turbine is a highly non-linear device, the Linearized 
state variable model developed in Ref. 2 is valid only for 
small changes from a specified steady state operating point. 
This model was obtained from a coupled ship and propulsion 
plant dynamic non-linear model consisting of the following 


features: 


A non-linear model of the LM2500 engine developed by 
General Electric. 

A non-linear model of the ship, propeller and propul- 
sion train provided by tne Bath Iron Works Company. 
Secondary flow losses considered. 

Bleed flow for extraction and turbine cooling considered. 
Inlet and exhaust losses incorporated. 

Pressure and temperature dynamics neglected. 

Nested loop balancing utilized to obtain various steady 


state operating points. 


In matrix equation form, this linearized state variable 
model of the various engine and ship parameters is expressed 
as: 
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CX + DU 


K 
ll 


X is the system state vector, and consists of three 


components: 
x, = NG - NG... ; NG = Gas Generator Speed (rpm) 
x, = NPT - NPT s ; NPT = Power Turbine Speed (rpm) 
X2 = V = Vos ; V - Ship Speed (knots) 


U is the system input vector, and consists of four 


components: 


u, = WES ES ; WF - Fuel Flow (lbm/hr) 

u, = β- Bos ; B = Stator Guide Vane Angle (deg) 

u, = WB - WB... ; WB = Compressor Bleed Flow (lbm/s) 
Dis P/D - PD. > P/D = Propeller Pitch/Diam Ratio ('"") 


Y is the system output vector, and consists of six 


components: 
Yan SMC = SMC -s ; SMC = Compressor Surge Margin (3) 
Og Τα = A ; Te = Combustor Outlet Temp (OR) 
Yan OPT - OPT, ; OPT = Power Turbine Torque (£t-1bf) 
Ya Pi P3 -s ; P3 = Compressor Total Disch Press (psia) 
Ys = HPPT - HPPT__; HPPT = Power Turbine Power (hp) 
Ye = T5 - De ; T5 - Power Turbine Inlet Temp en) 


and the ss subscripts indicate values at a specific steady state 
operating point, the same point used for determination of the 
constant A, B, C, D matrix elements. 

This study considered single engine up-power transients of 
70% to 75%, and 80% to 85% of maximum gas generator speed. In 
both cases, the steady state operating point selected was that 


REL 





for the higher power, or final condition. The necessary 
state variable data to accomplish this study was obtained 
trom computer printouts generated in Conjunction with Ref. 2, 


and is contained in Appendix B. Physical constraints on 


inputs and input rates are also contained in Appendix B. 
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III. CONTROL SYSTEM DSSIGN WITICCONSZN 


A. BASIC PHILOSOPHY 


The CONSYN (Control Synthesis) design program presented 
in Ref. 1 is a further refinement of linear regulator theory 
developed in Ref. 4. It is based on the initial assumption 
that the plant dynamics can be represented by the linearized 


state variable matrix equations 


PS 
li 


AX + BU 


[4 
| 


= CX + DU 


where, as previously discussed, all states are perturbed 

from their steady state operating points. This is illus- 
erated in LOCK diagram corm in Fig. las TRIS Original plant 
is then augmented with integrators to achieve open loop inte- 
gral control as shown in Fig. lb. Application of linear 
quadratic regulator theory to this augmented system leads 

to the minimization ot da performance or cost index, J, 


defined by the matrix equation 


Es 
, 


= (v-ov + ÜtmÜ) dt 


O 

where Q is a constant, diagonal, positive semi-definite 
matrix used to weight output deviation from final values, and 
R is a constant, diagonal, positive definite matrix used to 


weight input rate deviations from final values. 


I3 





The linear regulator problem solution yields an optimal 
state feedback regulator which follows the matrix equation 


control law 


Oy eG οἩ 


1 2 


and is illustrated in Fig. lc, where Z is the demand vector. 
Finally, in order to achieve integral control commands gen- 
erated by the error between demand and output, the system 
shown in Fig: lo is Converted tó that shown in Fig- ld. τις 
final system configuration has the same time response as that 
shown in Fig. lc but has the added advantage of assuring zero 
steady state error even if system degradation alters th 
elements of coefficient matrices. 

Since that choice of Q and R matrix elements for the 
solution of the linear regulator problem is arbitrary, 
COPES/CONMIN logic is then utilized to find the optimum Q and 
R which minimizes the cost index J without violating any of the 
imposed physical constraints. 

In summary then, the CONSYN design program applies a 
linear regulator design form to the given output control 


problem using the augmented state vector 
X 
Ke = [-1 . 


Block diagram methods along with suitably partitioned matrices 
are utilized to arrive at the desired structure of the optimal 


control system. 


14 





B. PROGRAM FLOW LOGIC 


The CONSYN design program operates in three distinct 
segments: Initial Design, Optimization, and Final Design 
Output. A summary of the significant logic and operations 


within each of these segments follows. 


Ly οι αἱ Design 


a. Input the following: 
([L) —A. B. C, D matrices 
(2) Initial guess for appropriate Q and R matrices 
( 3) Initial condition X and U vectors, which are 
then combined to form the augmented state 


vector K, = [-5-] 


(4) System constraints 
b. Calculate system open loop characteristics 


c. Calculate P, the solution to the steady state 


matrix Riccati equation. 


AP + PA, e pBR "Bep rco e. - 6 


| | 
where A = btg] B =C Cs =EjD] 
| | 


d. Calculate G the optimal state feedback gain matrix 


| 
Bi Y 
G=R B, P= LG) 1G] 


1 the detailed development of this Riccati equation and 
the following G, H, and L gain relationships in the partitioned 
matrix form is found in Refs. 1 and 4. 
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e, Calculate H and L the integral control gain matrices 


| -1 
[LH] = > η 


| 
where E = [1:10] and is required so that the number 
| 


of outputs fed back for comparison with the demand 
vector is equal to the number of inputs. 
f. Simulate time response of system. 


g. Print results of all of the above. 


2. Optimization 
Vary Q and R matrix elements under the direction of 
COPES/CONMIN logic to achieve the smallest cost function 
value without violating any constraints. 
3. Final Design Output 
Print complete results of system characteristics 
associated with the optimum Q and R matrix elements. 
Since the CONSYN program does not contain a plotting routine 
for system time response simulation, a separate program must be 
utilized. Forthis study, the IBM CSMP in conjunction with a 
VERSATEC plotter was employed. A derivation of the equations 


required for this simulation is contained in Appendix C. 


C. MODIFICATIONS TO THE ORIGINAL PROGRAM 


As discussed in Ref. 3, the following two significant modi- 


Fications were made to the original CONSYN program: 


The Kleinman technique for the solution of the matrix 
Riccati equation in subroutine OPGAIN was replaced by an 


elgenvalue method. 


16 





The fixed step integration solution in subroutine PEAK 
for system simulation was replaced by a variable step 


integration routine. 


In using CONSYN with these changes for this study, diffi- 
culty was encountered with the variable step integration routine. 
Numerical underflows, overflows, and step size decimation were 
encountered to varying degrees in all attempts to use it.  Con- 
sequently, the fixed step integration method was returned to 
subroutine PEAK, along with a logic modification to determine 
and print out the number of terms required for the matrix 
series convergence of PHI - S Here, the matrix F is 
defined as 


pon. dps 


where the matrices A,, B, and G are as previously defined. 
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Pig. las. Oraoumal Plant 





Piga Ib. Open Loop integral Control 


Figure 1. Integral Control Synthesis 
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Fig. le. Closed Loop Optimal Integral Control 





Fig. ld. Closed Loop Optimal Output Integral Control 


Figure l. (Cont'd) 


ES 





LV. RESULTS 


A. GENERAL 


The overall results were disappointing. Of all the 
basic design forms and variations evaluated, only one was 
found which did not violate one or more physical constraints. 
In addition, the study was hampered by the relatively long 
computer execution times required for the optimization of 
a Single proposed design (approximately 20 minutes), and the 
excessively long computer turn-around times (approximately 
10 hours). In all, at least 150 difference designs were 
attempted, with a total computer execution time expenditure 
of over 50 hours. A summary of the design forms and varia- 
tions considered is given in Tables I and II. 

Initial efforts to find an acceptable design were con- 
centrated on design form numbers 2 and 4 in Table I. These 
forms were selected because they appeared to be the most 
analogous to the design evaluated in Ref. 5 for the F401 
turbofan jet aircraft engine with respect to those outputs 
designated for feedback. All proposed designs Of these 
forms consisted of transients between 70% and 75% of maximum 
gas generator speed, variable bleed, and no state variable 


scaling. All variations considered resulted in violation 


See Appendix D for a more detailed discussion. 
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of one or more physical constraints. Examples most often 


noted were: 


Maximum propeller pitch exceeded 


9 


Maximum propeller pitch rate exceeded, either in the 


positive or negative direction 


Negative bleed 


Design form number 4 was then exercised at various 
transients between 80% and 85% of maximum gas generator 
speed with similar types of physical constraint violations. 

In an effort to determine why negative bleed was some- 
times occurring, G, J. Michael, one of the authors of Refs. 
4 and 5 was contacted. The following information and guidance 


was obtained from this communication: 


= 


Although variable bleed was used as a control in Ref. 


ἠ-- 


it also sometimes assumed negative values due to matne- 
matical anomalies in the engine model that could not be 
eliminated. As a result, bleed was not used as a 


control Inpub in Ref. >. 


Although the primary motivation for state variable 
scaling in Refs. 4 and 5 was to avoid publishing 
proprietary information in the open literature, doing 
so in the case of this study might eliminate additional 


computational anomalities. 
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TABLE 1 


DESIGN FORMS EVALUATED 


Basic Form I.D. Number 


__ E > ze EUN SO 
“a ua = uy si Ei 
U = 22 ο 32 23 B > 
E = 44 U4 u 2i 
"4 34 
u5 u- > 5η u5 ua 
u4 593 ^3 73 75 E 
Y3 Y3 76 Y6 Y6 *3 
: Y6 Y6 s *3 ο 21 
| Yi Ya Yi Yi Y, Y» 
Y > Y > Y» Y> Y> Y3 
Y4 Y, Y4 Y Y4 Ys 
Ys d fs 75 “3 76 


NOTE: The first n elements of the Y vector, where n is the 
number of inputs, are those that are to be fed back 
for comparison with the demand vector. 
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TABLE Il 
DESIGN FORM VARIATIONS EVALUATED 
Initial to final Power Level, % of Maximum Gas Generator 
Speed: 
a. 70 - 75% 


b. 80 - 85% 


a. Variable 

D. Constant 
Settling time - 20 to 200 seconds. 
Ship speed constraint with respect to settling time: 

a. Within * 29 of steady state value 

b. Not constrained 
State variable scaling or normalization with respect to 
the difference between initial and final steady state 
values: 

a. Scaled 

b. Unscaled 
Elements of the Q and R performance index weighting 
matrices allowed to be changed by the optimization 
routine: 

a. All Q and R elements 


5. ΤΊ (the weight of fuel flow rate error) fixed at 


1.0, all other Q and R elements allowed to change. 
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Design form numbers 5, 6, 7 and 8 were then exercised 
at various transients between 80% and 85% of maximum gas 
generator speed, constant bleed, and with state variable 


scaling.+ The following results were noted: 


An acceptable design which did not violate any of the 
input or input rate constraints was achieved with 
design form number 6, with the performance index 


weight on fuel flow rate error fixed at 1.0. 


There was no significant difference in system response 


obtained from design form numbers 5, 6 and 7. 


For those designs which did not impose a specification 
on ship speed with respect to requested settling time, 
the optimization routine reduced the degree of initially 
violated constraints during the first several intera- 
tions. The remaining iterations produced no furtner 
changes in the degree of constraint violation or in 


the value of the objective function. 


For those designs which did impose a specification on 
ship speed with respect to requested settling time, 
the degree of initially violated constraints was 
reduced at the expense of violating another constraint 


that was not initially violated. 


lsee Appendix E for the derivation of the scaled state 
variable data with constant bleed. 
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The number of terms required for the series convergence 
of PHI = ος in subroutine PEAK was much less than 
before state variable scaling was utilized; usually 


fewer than ten terms were needed for convergence. 


Deviations from actual values were noted for the 
linearized model calculations for outputs at the initial 
condition. The worst occurred for power turbine power 


which assumed a negative value. 


B. SPECIFIC 


A summary of the design data associated with the one 
acceptable design achieved is presented in Table III. Plots 
of system time response simulation for this design are shown 
in Figures 2a through 2f. Plots of state open loop response 
to a step input corresponding to the final steady state input 
values are shown in Figures 3a and 3b. It is noted that with 
the optimum design, all engine parameters reached their final 
steady state values within 6 seconds, and the ship speed 
achieved 99% of its final steady state value after 100 seconds. 

Summaries of performance data for various selected designs 
are presented in Table IV, with the degree of constraint 
violation indicated. 

A comparison of actual output values with those predicted 
by the linearized model at the initial condition is given in 


Table V. 
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TABLE ILE 


OPTIMAL CONTROLLER DATA, DESIGN 6-2F 


Performance Index Weighting Factors 


OULDUR, Error 





Ej 42 33 44 ds Ge ας dg 
(u,) (Y 4) (Y c) (x4) (yi) (y 5) (Y 4) (y 5) 
Noo .37  .26 7x10 ' 6x10 ας 72 .023 
The H Matrix 
13030 .5066 .9441 
1.533 ‚03452 = 026809 
797 9716 ος 
Mie L Matrix 
-.3585 -.02830 -.0006817 
.1184 2.095476 -.001542 
Bol so SS -.2850 


NOTE: 
this design is contained in Appendix F. 
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The complete computer output of the optimization run £ 


TnputgRateTeir:sr 


Be L > τι 
(à) (u) (u,) 
Tod) 32 28 





eonstraint 


ΜῈ 
max 


WF. 
min 


max 


min 
WB 
max 
WB . 
min 
P/D 


Dn 


NOTE: 


Limit 
7600 
900 


40 


© 


ὃ 
ὃ 


1413 
ο. 


+3200 
+59 


ho 


+.1 


Circled items indicate constraint limit exceeded. 


PERFORMANCE 





1229 


20:55 


-45.8 


ED 0 





34.7 
5) 
P 1 


TABLE IV 


DATA SUMMARY FOR VARLOUS DESIGNS 


1090 
10.2 


2x 








Design Code 
4=2B  6-2C 
3950 5970 
2360772360 
23.4 23.4 
1320051523 
190.939 ——--—- 
_ 

43 
12172212.05 
1920 GEL 

=12,7 -12,.7 


used above is explained as follows: 


(lst Digit) 


A 
i 


Design form 


from Table I 


= 


(2nd Dici) 
^ 
Power Code 


l sse» 
2 = 80-85% 
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6-2D 6=2E  6-2F 
4170 3860 3860 
2360 2360 2360 
241429 2559 7223.4 
ΕΠ 2 E 
constant = .02---- 
constant = .02---- 
ede lea eas 
a Το 
DO 1380 
eS €-17.2 ιο.” 


(letter) 


^ 
| 


VadrTeation Code 
from next pege 








8-2E 8-2Ε 
3950 3910 
2360 2360 
23.4 23.4 
porcum 
1.43 1.42 
E τε. 
3180 
CE m. 





TABLE: LV “(Cont da) 


Variation Code (Also see Table II) 


[> 


A eee 


Bleed 
variable x x 
Constant X X x x 


State Variables 
Scaled x x x 
Unscaled x X X 


Q and R Weights 
Variable x X X X X 
rj Fixed X 


Ship Speed 


Constrained X x x 
Unconstrained X x x 
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TABLE V 


COMPARISON OFS INITIAL CONDITION ACTUAL OUTEUTS WITH THOSE 
COMPUTED BY THE LINEAR MODEL 


Output Actual Value” Computed Value 
SMC 32659 30200 

T4 1776 Loos 

OPT 9706 8132 

P3 98.93 91.09 

HPPT 3273 -1795 

TS LLG 1365 

1 


Obtained from data generated in conjunction with Ref. 2. 


29 








E] € 
O 
Fae 
uJ 
m 
LJ 
j E 
C2 as 
z= =z 
O 
> 
5 
8 
O 
i 
| 
UN | 
O K 
ES a | 
7z | 
ἘΝ h 9 
E me 
am 
δα > 
δ = 
J mn 
4 l1 
C | = 
[1-1 | = 
ee 2 = 
«Ὁ \ 5 
= AS , 
CO 
PZN τ 
2) 
Dm 
— 5 
> τοὶ 
5) 
SOS 3 
Lf) im 
OO = 
| 
C=) i 
OS 2 
55 
00'0858  00'00h8 00°02E3 oo oneg 00-0918 ο ug 00°'0003 00'026824  90'Chg/7 
an 
—————————— 
OO He Omaha 00 0912 00 0002m 00 000e popa e CO’OMSI  00'09.1  QO'089! 
id 


30 


Figure 2a. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


CONCLUSIONS 


s 


Some form of state variable scaling is necessary to 
improve the efficiency of the numerical calculations. 
Scaling causes the A, B, C, D matrix elements to be 
much closer to the same order of magnitude. 

If bleed is to be utilized as a control variable, 

the failure of the mathematical model to prevent it 
from assuming negative values must be corrected. 
Although the system's non-linearity was in evidence 
at the initial condition state as Shown in Table Y, 
the mathematical model became a much more accurate 
representation of the actual system as the final 
steady state operating point was approached. 

Optimal Controller design with the CONSYN program is 
by no means a completely automated (or cooxbook) 
process. Much manual interaction with options and 
variable parameters must be employed in order to 
achieve a feasible design. 

Input and input rate constraints were design limit- 
ing. As long as these constraints were not violated, 
there were no unreasonable excursions of any state or 


output variables. 
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B. RECOMMENDATIONS 


l. The specific reason(s) why the optimization routine 
is so frequently unable to reduce violated constraints 
beyond a certain point must be investigated and 
corrected before any practical controller design by 


this method can be achieved. 
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Name 


ANALIZ 


CONTL 


COST 


DET 
INVERT 


LIAP 


MULT 
MULTV 
MULTZ 
OBSERV 
OPGAIN 
PEAK 
PHIOFT 
POLY 
PRINTM 
PRINTN 
RANK 


ROOT 


APPENDIX A 


LIST OF CONSYN SUBPROGRAMS 


Purpose 


Integral control design and analysis (pre- 
viously named CYCLE in Ref. 1) 


Controllability Check 


Evaluation of Performance index and deter- 
mination of maximum eigenvalue 


Matrix determinant 
Matrix inverse 


Solution to bilinear symmetric Lyapunov 
equations 


Matrix multiplication 

Vector multiplication 

Complex matrix multiplication 
Observability check 

Solution to optimal linear regulator 
System simulation 

State transition matrix 
Coefficients of characteristic polynomial 
Matrix print, double precision 
Matrix print, single precision 
Matrix rank 


Polynomial roots 
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1. State Variable Representation at 703, 


APPENDIX B 


PROPULSION PLANT DATA 


79$ 


80%, and 85% of Maximum Gas Generator Speed 


States 


NG (rpm) 
NPT (rpm) 
V (knots) 


Inputs 


WF (lbm/hr) 
8 (deg) 

WB (1bm/S) 
η) 


Outputs 


SMC (8) 

T4 (OR) 

QPT (ft-lbf) 
P3 (PSIA) 
HPPT (hp) 

T5 ("R) 
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APPENDIX B (Cont'd) 


2. A, B, C, D MATRIX ELEMENTS FOR 75% 
OF MAXIMUM GAS GENERATOR SPEED 


The A Matrix 


=. 21709 =+1186 ου 
«1295 -.9364 Sd ial 
Qi... .0002429 zug 3 


The B Matrix 


1281 37.74 = A 0.30 
08926 “το πο ο -328,6 
00 ο ου -1114 


The C Matrix 


-01073 =.000327L 0.0 
-.2160 «00081984 9 
SITE -2.473 00 
ο... 2000 1613 0:0 
1928 4526 0.0 
-.2289 .01866 0.0 
The D Matrix 
-.01811 -.5678 2.453 00 
.4905 9.434 τα”. 0.30 
25120 -219.9 -192.5 0.0 
.008820 -1.405 =L.L95 050 
33.00 -54.97 -48.13 QU) 
3035 11.49 11.85 0.0 
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3: A, B3 


The A Matrix 


= bee 
. 3466 
oTo 


The B Matrix 
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The C Matrix 
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The D Matrix 
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APPENDIX B 


(Conte a) 


C, D MATRIX ELEMENTS FOR 85% 
OF MAXIMUM GAS GENERATOR SPEED 
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APPENDIX B (Cont a) 


4, SCALED A, B, C, D MATRIX ELEMENTS FOR 85$ 
MAXIMUM GAS GENERATOR SPEED, CONSTANT BLEED 


The Scaled A Matrix 


2212 u.s 
«905 ΙΙ 
0-0 02029 


The Scaled B Matrix 
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The Scaled C Matrix 
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APPENDIX B (Cont'd) 


PHYSICAL CONSTRAINTS 


Note 
Inputs 
900 < WF < 7600 lbm/hr (1) 
SO (1) 
OL < P/ D < I3 (2) 
Input Rates 
-3200 « WF « 3200 lbm/hr/S (3) 
EEE: Er 
SD DASS (4) 
NOTES: (1) Approx. range between idle and full power 
Condition: 
(2) Model does not allow for negative pitch. 
(3) Obtained TOn C.E; 
(4) Computed from data provided by Bird Johnson 


Company as follows: 


Propel ler Diam er = LO. See 
Maximum Ahead Pitch =. 23.45 ει 
Maximum Reverse Pitch = 14.7 ft 


Minimum Time for Pitch Change from 
Full Reverse to Full Ahead = 20 sec 


A P/D _ 23.5 _ | 


AE 153 
20 5 


Pu 


er, = 2 
Τετ) ο S 





P/D = 
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APPENDIA € 


DERIVATION OF EQUATIONS USED 
FOR SYSTEM SIMUBATION WITH CSMP 


The following matrix equations may be obtained from 


inspection of the system block diagram in Fig. ld: 


X = AX + BU 


K 
li 


C 
il 


O 
U = 


However, 


desired output deviations are zero at the final steady 


state op 
Omg in: 
Ú = 
Recallin 
E = 
the abov 
into the 
g-output 


zo s 


x, = 


EX = ου 


S[H(Z-EY) Jät - LX 

Fr, UPON GlETerentiatand, 

H(Z-EY) - LX 

the demand vector Z is equal to zero since all 

erating point which has been defined as the 
Thus, 

- (HEY + LX) 

q that the matrix E was defined as 

[110], 

e matrix equations can, therefore, be expanded 
following form for the 3-state, 3-input, and 
system considered: 


a Xy T 812 x + a 


I 2 l3 


X +a 


57. 33 


12 τ 


81 ο ο ο δα. 
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APPENDIX D 


RATIONALE FOR DESIGN FORM SELECTION 


In Ref. 5, the following outputs were selected for 
feedback in the case of the F401 turbofan jet engine 
controller design: 

l. Fan inlet guide vane angle (FIGV) 

2. Rear compressor variable vane angle (RCVV) 

3ed Κατι 

4. High turbine inlet temperature (HTIT) 
Those outputs selected for feedback in design form number 
2 of this study were: 


l. Propeller pitch to diameter ratio (?/D) 


2. Compressor stator guide vane angle (3) 
3. Power turbine torque (QPT) 
4. Power turbine inlet temperature (T5) 


A comparison of these designs follows: 

1. The marine propulsion gas turbine engine, of 
course, has no turbofan. However, the P/D was 
considered analogous to the FIGV since both 
directly relate to the loading of the output 
turbine. In the aircraft case, the output 
turbine drives the fan; in the marine case, 
the output turbine drives the propeller. 

2. B is completely analogous to RCVV since this is 


the compressor vane control in both cases. 
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APPENDIX D (Cont'd) 


The thrust in the aircraft case 1s a measure of 
the total power produced. In the marine case, 
either ship speed, power turbine power, or power 
turbine torque is analogous. Torque was selectec 
since it was considered the easiest of the three 
to measure accurately. 

HTIT is equivalent to combustor outlet temperature. 
However, the nigh temperature involved shortens 
transducer lifetime, and is variable around the 
perimeter of the combustor. To eliminate tnese 
difficulties, power turbine inlet temperature, 


instead, was chosen for feedback. 
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APPENDIX E 


DERIVATION OF SCALED STATE VARIABLE DATA 
WITH CONSTANT BLEED 


From the matrix equation, 
X = AX + BU 
consider one of the scalar equations 


TAD 


tara * 844 X4 * 2941; U, * 515 U5 


ΤΠ 


+ D ur. + b 


¡cs pq ed 


t 
[] S . « 
Since the x and u = are deviatıons from the steady 


state Operating point, and bleed 1s constant, then u. = Ox 
Now define the following scaled variables: 


τ 1 
X. ==> ἂ. = 


where the AS refer to the difference between initial and 
final conditions, and Uy will not be scaled since its 
initial and final values are equal. 

Then 


» 


ra AX, X, E 81 3 2X3 X3 


AX- -AX 12 


"Τι en, 
+ 2411 Au, uy ~ D1 5 LU, U, - 214 Uy 


Dividing through by AX, results in 





a AX, , aa Au, | 
ΟΥ, sen 3 πε p 
Au D 
ο. 14 - 
+ b, — u, + — 
12 AX, 2 âX] 4 
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APPENDIX E (Cont'd) 


Carrying this procedure out with the other state and 


output equations results in the scaled matrix equations: 


^ 7 NN 


x = AX + BU 

















Y = CX + DU 
where 
Ν i AX. i 2X4 
2 E = E ux 
TL 1 Xy 13 ax, 
A = a 823 8 — 
21 Ax) 2 23 AX, 
] AX} AX. 
30M. — 73999 — 33 
3 3 
2 Au, : Au, 4 1 
A i A: 
TUE AX, 12 AX, 14 AX, 
B = b σπα b i CER 
2L AX, 22 AX, 24 AX, 
7 Au, E Aus 4 1 
SL Ax 22 AX, 34 AX, 
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COMPUTER OUTPUT OF OPTIMEZATION 
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